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ABSTRACT 

Some  nonlinear  stationary  reaction-diffusion  systems  involving  nonlinear 
terms  which  may  be  discontinuous  are  considered.  Such  systems  occur,  for 
instance,  in  the  study  of  chemical  reactions  and  the  discontinuities 
correspond  to  reactions  of  order  zero.  In  such  concrete  model,  the  set  where 
the  reactant  vanishes  plays  an  important  role.  Here  we  prove  the  existence  of 
solutions  for  a  general  class  of  such  systems  satisfying  Dirichlet  or 
nonlinear  boundary  conditions.  Necessary  and  sufficient  conditions  are  given 
assuring  that  the  reactant  component  vanishes  on  a  set  of  positive  measure. 
Estimates  on  the  location  of  such  set  are  given. 
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SIGNIFICANCE  AND  EXPLANATION 


Nonlinear  reaction-diffusion  systems  arise  in  a  variety  of  applications 
including  mathematical  biology,  chemical  reactions,  and  combustion.  The 
nonlinear  terms  in  the  model  equations  are  often  merely  continuous  or  even 
discontinuous.  A  typical  example  occurs  in  the  study  of  a  single  irreversible 
nonisothermic  reaction  of  order  less  than  one  and  the  discontinuous  terms 
correspond  to  reactions  of  order  zero. 

In  this  paper  the  authors  formulate  the  model  nonlinear  system  with  the 
help  of  maximal  monotone  graphs  and  they  prove  the  existence  nf  solutions  for 
a  general  class  of  stationary  systems  (including  the  model  of  chemical 
reactions)  satisfying  Dirichlet  or  nonlinear  boundary  conditions.  The 
regularity  and  uniqueness  of  solutions  are  also  considered.  Of  particular 
interest  is  the  existence  (or  nonexistence)  of  a  dead  core  i.e.  a  set  in 

which  the  reactant  component  vanishes.  This  set  plays  an  important  role  in 
the  modelling  of  chemical  reactions.  Necessary  and  sufficient  conditions  are 
given  for  the  existence  of  as  well  as  estimates  for  the  location  of  such 


a  dead  core  whenever  it  exists. 
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J.  Ildeforiso  and  Jesus  Hernandez ^ 

40.  Introduction. 

3  3  3333SauM 

Many  papers  have  been  devoted  during  recent  years  to  the  study  of 
reaction-diffusion  systems  which  arise  very  often  in  applications  such  as. 
e.g. ,  mathematical  biology,  chemical  reactions,  combustion  theory,  etc. 

Here  we  consider  a  system  describing  a  single,  irreversible, 
nonisothermic  stationary  reaction  of  the  form 

Y-I 

f  -Au  +  y2F(u)e  v  *  0  in  Q 


(0.1) 


Y-I 

2  '  v 

-Av  -  vu  F(u)e 


0  in  Q 


du 

*5 —  +  €  (u-1  )  »  0 


on  32 


(0.2) 


+  C(v-1 )  «  0 


on  32 

where  2  is  a  bounded  open  subset  of  E1*,  ji2  is  the  Thiele  number,  v  is 
the  Prater  temperature  and  Y  is  the  Arrhenius  number  (see  [3]).  Here  e 
and  i  (the  Biot  numbers)  are  positive,  being  in  some  cases  infinity,  in 
which  case  (0.2)  is  interpreted  as  the  Dlrichlet  boundary  conditions 
(0.3)  u  =  1  ,  v  ■  1  on  32. 

The  function  F(u)  is  assumed  to  be  nondecreasing  and  it  is  also  assumed  to 
satisfy  F(0)  «*  0,  F(  1 )  *  1  and  F(s)  >  0  if  s  >  0.  The  unknowns  u  and 
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v  are  nonnegative  and  they  represent,  respectively,  the  concentration  and  the 
temperature  of  the  reactant. 

Very  often  F  takes  the  simple  fora  F(u)  -  vP,  where  p  >  0  is  the 
order  of  the  reaction  (see  [3]  Vol  I).  In  the  case  of  a  reaction  of  order 
zero  F  is  given  by  F(0)  -  0  and  F(s)  *  1  if  s  >  0  (thus  F  is  a 
discontinuous  function). 

Existence  and  uniqueness  results  for  the  parabolic  problem  associated 
with  (0.1),  (0.2)  (or  (0.1),  (0.3))  have  been  given  by  some  authors  (c.f.  e.g. 
[2],  [4],  [20]).  Existence  and,  in  some  particular  situations,  uniqueness 
results  for  the  elliptic  problem  can  be  found  in  [2],  [23],  [17]  and  [21]  for 
p  >  1.  The  case  0  <  p  <  1  is  considered  in  [3]  p.  311  (see  also  [22])  but 
existence  theorems  are  not  given.  It  is  shown  in  [3]  and  [22]  that  for 
p  ■  0,  if  m  is  large  enough,  no  strictly  positive  solution  can  exist.  It 
is  also  shown  that,  in  some  particular  examples,  the  set 

■  {x  €  0  s  u(x)  «  0}  (called  the  dead  core)  is  not  empty  and  has  positive 
measure  if  0  <  p  <  1. 

The  main  idea  used  in  [22]  and  many  other  papers  (cf.  [3])  is  to  reduce 
(0.1),  (0.2)  to  a  nonlinear  elliptic  boundary  value  problem  for  u  alone. 

Here  we  follow  a  different  approach  which  allows  us  to  obtain  better 
results.  Moreover,  we  are  able  to  treat  the  case  of  nonlinear  boundary 
conditions,  which  cannot  be  handled  by  the  preceding  device. 

We  shall  consider  the  case  of  discontinuous  functions  F(u)  in  the 

framework  of  maximal  monotone  graphs  in  ft2  (see  [8]).  For  the  reader's 

2 

convenience,  we  recall  that  a  maximal  monotone  graph  a  in  K  is  always 
specified  by  a  real  non-decreasing  function  6  by  a(r)  -  (-*»,8(r)]  if 
8<r)  -  a(r)  -  [6 ( r- ) , 8 ( r+) ]  if  -•  <  3(r-)  <  8(r+)  <  +-  and 
o(r)  -  [8(r-),4»)  if  8(r+)  -  +».  We  define  D(o)  -  {r  «  R  :  a(r)|«  0}  and 


a 
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the  sections  a  and  a  by 


a+(r)  -  max{z  s  z  e  a(r)}  if  re  D(a)  , 

a~(r)  -  mintz  s  z  9  a(r)}  if  re  D(a)  , 

a+(r)  =*  a”(r)  ■  ¥•  if  r  $  D(a),  r  >  sup  D(a)  , 

a+(r)  -  a‘(r)  -  --  if  r  $  D(a),  r  <  inf  D(o)  . 

Finally,  we  define  a°(r)  as  the  element  of  a(r)  with  minimal  absolute 
value. 

Through  the  paper  we  shall  study  the  following  genera-l  formulation 
including  the  system  (0.1)  as  a  particular  case: 


(NLS) 


f  -Au  +  a(u)f(v)  9  0 
V  -Av  -  0(u)g(v)  9  0 


in  0 
in  0 


with  the  boundary  conditions 

(DBC)  u  »  ^  ,  v  »  on  M 

as  well  as  the  nonlinear  boundary  conditions 


(NBC) 


_  3u 

s  TK 

_  3v 
=  ^ 


+  b(u) 
+  c(v) 


on  38 
on  38 


where  8  is  a  bounded  open  subset  of  B?1  with  smooth  boundary  38.  We  also 
assume  for  the  rest  of  the  paper  that 

(0.4)  a  and  0  are  maximal  monotone  graphs  such  that  0  e  a(0)  n  0(0)  . 

(0.5)  f  and  g  are  C1  functions  and  f(s)  >  0,  g(s)  >  0  if  s>0  . 

(0.6)  wt,  +,  and  t2  6  C2(3D) 

(0.7)  b  and  c  are  C  non-decreasing  real  functions. 

In  particular,  if  a  and  0  are  single-valued  (i.e.  they  are  continuous  real 
functions)  then  the  set  inclusion  of  (NLS)  should  be  replaced  by  equality. 


•mmmam 


nmtsM 


In  the  general  situation,  (u,v)  6  H2(ft)  *  H2(0)  is  a  solution  of  (NLS) 
if  there  exists  a,  d  6  L2(Q)  such  that  a(x)  6  a(u(x)),  d(x)  S  0(v(x))  a.e. 
x  S  Q  and 

-Au  +  af(u)  *  0  and  -Av  -  dg(v)  *  0  in  0  ■ 

Me  shall  prove  the  following  existence  result,  which  extends  in  some 
sense  those  in  (2],  (23]  and  (17). 

Theorem^A*  Assume 

(A.1 )  D(o)  -  D($)  -  R  , 

(A. 2)  f  ( S )  >  «1  >  0  VS  e  R  , 

and 

(A* 3)  0  <  g(s)  <  m2  Vs  6  R  . 

Then  there  exists  at  least  one  solution  (u,v)  of  (MLS)(DBC)  (resp. 

(MLS)  (NBC) ).  Moreover  u,  v  6  W2'p(ft)  for  any  p,  1  <  p  <  +**. 

Me  also  consider  the  existence  and  non-existence  of  a  dead  core  ftQ 
where  u  *  0  and  consequently  the  existence  of  the  free  boundary  ^ . 

Roughly  speaking,  such  a  dead  core  for  (0.1),  (0.3)  arises  when  it  is 
Impossible  for  diffusion  to  supply  enough  reactant  from  outside  ft  to  reach 
the  central  part  of  ft.  (cf.  [22]).  This  may  happen  if  the  reaction  rate 

y-i 

F(u)e  remains  high  as  the  reactant  concentration  decreases.  Thus  (for 

(0.1),  (0.3))  the  existence  of  ftQ  depends  essentially  on  three  things:  the 
reaction  order,  the  Thiele  number  and  the  size  of  ft. 

Our  main  result  in  this  direction  can  be  stated  in  the  following  general 
terms. 


(1) 

There  is  a  large  literature  about  this  subject  in  the  case  of  a  single 
nonlinear  equation.  See,  e.g.  the  systematic  study  of  [12]. 
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Theorem  B.  Assume  that  die  hypotheses  of  Theorem  A  are  satisfied*  Then  the 
following  properties  are  true: 

i)  Tf_  <*(s)  *  u2|s|p"1s  and  (u,v)  is  any  solution  of  (NLS)(DBC)> 
then  a  dead  core  may  exist  only  if  0  <  p  <  1.^ 
ii)  Let  o(s)  *  W2|s|P_13  with  0  <  p  <  1  and  let  (u,v)  a  solution 
of  (NLS) (DBC)«  For  X  >  0  let 

=  tx  6  n  :  f(v(x))  *  x} 


(0.8) 


Then 

IlR 

3  {x  S  Q  }  d(x,3Q.  -  OQ  -  supp  *  ))  >  [— — )  2  } 
0  \  X  1  KX,U 

where  M  =  l^.l  _  and 

1  L  00) 


iii)  Let  o(s) 
M  b£  M 


K 


X,y 


li2*sign  s. 

Izl  „  , 

l  (n> 


i 

f2(1+p+(N-1  )(p-1  )]  'P-1 
‘  v  2  2  '  * 

Xp  (1-p) 

Then  the  estimate  (0.8)  holds  if  we  replace 

2 

where  z  satisfies  Ai  *  P  n1  in 


Q  and  z  ■  ^  on  3Q. 

Furthermore,  if  n  is  convex,  the  above  results  are  still  valid  for 

(NLS) (NBC)  in  the  sense  that  if  0  <  p  <  1 ,  then  n^  has  a  positive  measure 

for  |i  large  and  it  is  possible  to  estimate  QQ  (see  (2.22)). 

The  above  theorem  is  specially  meaningful  if  m.,  in  (A. 2)  is  strictly 
positive  (this  is  true  in  the  case  of  the  combustion  system  (0. D),  (0.3)) 
because  then  we  have  0^  ■  n  for  any  X  8  (0,m^ )  and  so  the  estimate  (0.8) 
reads 

Air 

(0.9)  Q  3  (x  8  Q  :  d(x,3n)  >  2  }  . 


By  convention  |s|”1s  ■  sign  s  (*  -1  if  a  <  0,  *  if  s  *  0  and 

-1  if  s  >  0). 
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From  the  definition  of  K. 


in  Theorem  B  we  deduce  that  K,  ^  0  when 

A  j|J 

X  ^  0  or  v  0,  and  that  K,  *  +•  if  X  P  +»  or  n  ^  Therefore  for 

1 

a  fixed  bounded  Q  the  existence  of  a  dead  core  nay  only  be  guaranteed 

(by  estimate  (0.9))  if 

±-E 

diameter  of  Q  >  (•=-^ — )  2  , 

m,,)! 

assuming  0  <  p  <  1.  Then  a  critical  value  |i  of  u  can  be  found  such  that 

c 

£1  has  positive  measure  if  u  >  u  .  In  fact,  direct  computations  show 
0  c 

(when  N  -  1 )  that  function  u  is  strictly  positive  in  (l  if  u  <  u  , 

■Bsassaasas  C 

and  u  vanishes  only  at  one  point  if  u  ■  |i  (see  the  proof  of  Lemma  2.1  and 

c 

also  [22]).  Estimate  (0.8)  of  Theorem  B  can  be  also  written  independently  of 
the  function  v  for  other  systems  in  which  it  is  not  difficult  to  estimate 
the  set  (for  instance  *  Q  if  in  (MLS)  we  assume  f(s)  =  s  and 

2  >  6  >  0  for  $  large  (or  X  small)  enough). 

Through  the  paper  we  also  remark  on  other  more  general  formulations  of 
(NLS).  The  parabolic  problem  associated  with  (NLS)  will  be  studied  in  a 
forthcoming  paper  by  the  authors  [13].  The  case  of  0  unbounded  will  be  also 
treated  elsewhere  [14]. 
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§1.  Existence  results. 

Consider  first  the  problem 


(DP) 


-Au  +  a(u)f(v)  a  0  in  0 
-Av  -  fHu)g(v)  8  0  in  Q 
u  ■  ^  ,  u  *  V  2  on  30  , 


where  0  is  a  bounded  open  subset  of  with  smooth  boundary  30  and 

a,  P#  £»  g ,  ,  <P 2  satisfy  (0.4),  (0.5)  and  (0.6).  Set  X  -  (H2(0))2. 

Definition  1.  We  shall  say  that  (u,v)  6  X  is  a  solution  of  (DP)  if  there 

nflaaaaassan  ■  ■  . .  1  **  . . . 

exist  functions  a,  b  0  L2(fl)  such  that  a(x)  6  a(u(x)),  b(x)  6  0(v(x))  e.e. 
in  0  and 

-Au(x)  +  a(x)*f(v(x))  -  0  a.e.  x  8  0  , 

-Av(x)  -  b(x)g(v(x))  •  o  a.e.  x  e  0  , 


and  the  boundary  conditions  (DBC)  are  satisfied. 


Definition  2.  The  pair 

((un,vn),(u°,v0)]  8  x  x  x 

is  a  sub-supersolution  < 

naaaaa: 

(DP)  if 

*  o  *  0 

.  uo  *  u  *  vo  *  v 

a.e.  on  0  and 

(1.1) 

-Auq  +  a“(uQ)f(v) 

4  0  4  -AuQ  +  «+(u°)f(  ) 

W  C  [v0,v°]  , 

(1.2) 

-Avq  -  B°(u)g(vQ) 

4  0  4  -Av°  -  B°(u)g(v°) 

Wi  8  lu0,u°l  , 

(1.3) 

uo  <  *i  *  u° 

on  30  , 

(1.4) 

vQ  <  V»2  <  V° 

on  30  , 

where 

tK,i)  »{h8  l2(0)|k(x)  4  h(x)  4  t(x)  a.e. 

_on  0}  if  K,i  8  D2(Q). 

Our  main  existence  result  for  (DP)  is  the  following 
«22E2?«iI*  Suppose  that  [  (uQ,  vQ),  (u°,v°)  1  is  a  sub-sni»ec3olution 
satisfying 

(H,)  u0,v0,u0,v°  8  d"(0) 


and  that 


D(a)  =  D(&)  »  R 


(h2) 

Then  there  exists  at  least  one  solufon  (u,v>  of_  (DP)  such  that 
uQ  <  u  <  u°,  vQ  <  v  <  v°.  In  addition  u,v  6  W2'^((l)  for  any  p# 

1  <  p  <  4». 

Remark  1.1.  This  theorem  generalizes  results  of  [2],  '17],  (18)  and  (22). 

2  2 

To  prove  Theorem  1.1  we  define  E  *  [L  (Q)]  and 
K  =  [uQ,u°]  *  [vQ,v°].  It  is  clear  that  K  is  a  convex,  closed  and  bounded 
subset  of  E.  Now  we  define  a  nonlinear  operator  TsK  ♦  E  in  the  following 
way:  for  (u,v)  6  K,  T(u,v)  -  (w,z)  is  the  unique  solution  of  the  uncoupled 
system 


(1.5) 

-Aw  +  a(w)f(v)  +  w  9  u 

in 

n 

(1.6) 

w  - 

on 

an 

(1.7) 

-Az  +  M.z  *  8®(u)g(v)  + 

M.V 

in 

n 

(1.8) 

w  -  V>2 

on 

an 

• 

Here  M 

>  0  is  such  that  the  right-hand 

side 

of 

(1.7) 

is  increasing  in  v 

(we  can 

choose  such  a  M  because  g  is 

C1 

and 

(Hn) 

has  been  assumed). 

Indeed  by  (H2)  we  can  apply  the  results  of  Brezis-Strauss  [10]  to  obtain  the 
existence  of  a  unique  solution  w  of  [1.5],  (1.6).  Moreover,  by  (H^ ) ,  (H2) 
and  the  Lp-regularity  results  (see  e.g.  [16])  w  6  W2'p(fl)  for  any  p, 

1  <  p  <  •♦*».  A  similar  argument  works  for  z. 

The  proof  of  Theorem  1.1  will  follow  from  Schauder's  fixed  point  theorem 
applied  to  the  operator  T.  It  is  sufficient  to  check  that  T  is  compact  and 
that  it  sends  K  into  itself. 

Lemma  1.1.  T  is  compact. 

>394*3393  ■  ■*■■■■ 

Proof .  As  K  is  bounded  it  is  easy  to  show  that 

«wl  <  C 

with  C  independent  of  (u,v)  6  X.  Thus  it  is  suf ficieui-.  to  recall  the 
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1  2 

compactness  of  the  imbedding  H  ((2)  ^  L  ((2)  to  see  that  T  sends  bounded 
subsets  into  relatively  compact  ones  (the  same  for  z).  To  prove  that  T  is 
continuous,  suppose  that  ^un'vn^  *  (u,v)  in  E.  Then 

-A(w-w  )  +  o(w)f(v)  -  a(w  )f(v  )  +  w  -  w  3  u  -  u  in  (2 
n  n  n  n  n 

w  —  w  =  0  on  9(2  . 

n 

Multiplying  by  w  -  wn  and  integrating  by  parts  we  obtain  (for  the  case  a 
singlevalued  for  simplicity) 

/n|V(w-w  )|2  +  /0[a(w)f(v)  -  o(w  )f(v  )l(w-w  )  +  /0|w-w  |2  » 

■  fn I V (w-w  )|2  +  /_[a(w)f(v)  -  a(w  )f(v))(w-w  )  + 
u 1  n  si  n  n 

+  fa  o(wn>(f(v)  -  f(vn))(w-wn>  +  Jfllw-wn|2  *  /jj(u-un)  (w_wn^  > 

>  /fl|V(w-wn)|2  +  /n  a(wn)(f(v)  -  f(vn))(w-wn>  , 

and  by  the  Cauchy-Schwart  inequality  it  follows  that 

Iw-w  I2,  <  la(w  )l  _  If ( v)  -  f ( v  )l  ,  lw-w  I  _  + 

nH1((2)  "  L  (8)  "  L2((2)  n  L2((2) 

+  lu  -ul  Iw-w  I  . 

n  L2((2)  n  LZ((2) 

Now  it  is  easy  to  conclude  that  w^  ♦  w  in  H^{8).  A  similar  argument  can  be 
used  for  z. 

Lemma  1.3.  T(K)  c  K 

S3 3333333 

Proof .  We  first  prove  uQ  <  w,  i.e.  (uQ  -  w)+  =  0,  with  h+  =*  max(h,0). 

For  v  *  v,  (1.1)  yields 

0  >  -A(Uq-w)  +  o(uQ)f(v)  -  u(w)f(v)  +  uQ  -  w  . 

(We  again  suppose  a  singlevalued  for  simplicity  in  the  notation.)  Multiply 
this  inequality  by  (uQ  -  w)+,  integrate  over  (2  and  use  Green's  formula  to 
obtain 
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0  >  fa  -A(u0-w)(uq-w)+  +  fa  (a(uQ)  -  a(w))f(v)(uQ-w)+  + 

♦  Iq  <“0-w)(,Vw)+  *  /qIV(uo"w)+|2 

by  the  aonotonicity  of  a.  This  gives  (uQ-w)+  -  0.  A  similar  argument  shows 
that  w  <  u°. 

For  the  second  component  v  we  have,  with  u  **  u  in  (1.2), 

0  >  -A(vQ-z)  +  3°(u)g(v)  -  3°{u)g(v°)  +  M(vQ-z)  . 

Multiplying  by  (vQ-z)+  and  integrating  yields 

0  >  -fa  A(vQ-z) (Vq-z)+  +  /nt$°(u)g(v0)  +  MvQ  -  B°(u)g(z)  -  Mz)(v0-z)+  > 

(by  the  choice  of  M  the  second  integral  is  positive). 

Then  T  has  at  least  one  fixed  point  (u,v)  in  K  which  is  a  solution 
of  (DP).  Moreover  u,  v  6  L  (ft)  and  this  implies  that  u,  v  6  wi,p(ft)  for 
any  p,  1  <  p  <  **. 

Remark  1.2.  It  follows  easily  from  Morrey’s  theorem  (w2,p(ft)  ***-c1,r(ft) 

if  p  >  N  with  r  *  1  -  j^)  that  u,  v  8  C1r^(0)  for  any  0  <  6  <  1.  On  the 

other  hand,  if  we  suppose  for  instance  that  a  and  3  are  C1  then  u, 

2  6 

v  6  r'  (S)  for  every  0  <  6  <  1.  Indeed,  in  this  case  a(u)f(v), 

B(u)g(v)  8  C^(fi)  and  we  can  apply  Schauder  theory  ([16]). 

The  main  conclusion  of  Theorem  A  (for  the  Dirichlet  problem)  follows  from 
the  next  lemma. 

Lemma_1 .4.  Suppose  (H1 ) ,  (H2)  and 
(H3)  0  <  ra1  <  f ( 8 )  Vs  9  R  . 

Then  if  uQ,  vQ,  u°,  v°  8  H2(ft)  satisfy 

(1.9)  “Auo  +  mia""^uo^  <  0  *  “Au°  +  m1a+(u°)  ft 

(1.10)  Uq<^1<u°  on  3ft 

(1.11)  -Av  -  3°(-lw  I  )g(vn)  <  0  <  -Av°  -0°C»^t«  • 

0  1  L  (3ft)  0  1  L  (30) 
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)g(v°)  in  0 


(1.12) 


v0  <  *2  <  v 


on  3Q 


the  couple  ((Uq,Vq),  (u°,v0)]  is  a  sub-aupersolution  for  (DP). 

*  0  *  0 

Proof.  Let  viQ<u  <  u  ,  vQ  <  v  <v.  By  the  maximum  principle  we  have 

-  I*  I  _  <  m  <  0  <  u°  <  W  I 

1  L  (30)  °  1  L  OQ) 


Then,  by  (1.9) 


-Auq  +  o  (uQ)f(v  )  <  -Auq  +  n^a  (uQ)  <  0 
-Au°  +  a+(u°)f(v*)  >  -Au°  +  m1o+(u°)  >  0  , 
and  also  by  (1.11) 

-Av  -  B°(u*)g(v  )  <  -Av  -  I  )g(v  )  <  0 

°  °  1  L  OQ)  0 

-Av°  -  B°(u  )g(v°)  >  -Av°  -  B°(l^.l  m  )g(v°  >  0  . 

L  on) 

Moreover,  a  simple  argument  gives  v^  <  0  <  v°. 


0  0  • 

Now  the  problem  is  to  find  uQ,  u  ,  vQ,  v  6  L  ( n )  satisfying  (1.9)  - 

(1.12).  The  fact  that  such  uQ,  u°  and  vQ  exists  follows  from  the  results 

of  [10]  applied  to  a  and  B  respectively  (we  point  out  that 

I  )  <  0).  Hence  we  only  need  a  solution  v°  of 

L  on) 


-Av°  -  B°(l<P  I  }g(v°)  >  0  in  n 

L  Ott) 


v°  >  * 


on  an  . 


It  is  clear  that  it  is  sufficient  to  take  as  v°  the  (unique)  solution  of  the 


linear  problem 


-Aw  ■  ai#2  in  n 


w  ■  i on  3n  . 


This  proves  Theorem  A 


Remark  1.3.  It  is  clear  that  assumption  (A. 3)  of  Theorem  A  is  only  used  to 
find  v°.  if  for  example,  g  is  such  that  the  nonlinear  problem 

-Aw  ■=  8°(lv>  I  „  )g(w)  in  0 

1  L  (30) 

w  *  on  30 

has  a  solution,  then  we  can  remove  (A. 3).  There  is  a  very  extensive 
literature  for  this  kind  of  problem  with  different  assumptions  on  g,  but  we 
do  not  want  to  consider  this  point  here  (cf.  e.g.  [1]  and  the  survey  [19])* 
Remark  1.4.  If  a  is  assumed  to  be  singlevalued  and  C1  the  hypothesis 
f(s)  >  m2  >  0  is  not  necessary  (cf.  [17]).  On  the  other  hand,  if  a  and  0 
are  singlevalued  and  a,  0,  f  and  g  are  C1  with  sufficiently  "small" 
Lipschitz  constants,  then  it  can  be  shown  (cf.  [2],  [17])  that  the  solution  is 
unique. 

It  is  very  easy  now  to  apply  the  preceding  results  to  the  particular 
example  (0.1),  (0.3)  considered  at  the  beginning  of  this  paper.  It  is 

y-  — 

sufficient  to  take  f(v)  ■  g(v)  ■  e  v,  o(u)  -  M^up,  0(u)  ■  vu2up,  p  >  0 
and  *  *2  m  1  *  A  sub-supersolution  is  given  by  uQ  =  0,  u°  =  1 ,  vQ  =  0 
and  v°  the  unique  solution  of 

0  2  T 

-Av  ■  vy  e  in  0 
v®  ■  1  on  30  • 

The  case  of  nonlinear  boundary  conditions  can  be  handled  in  a  very 
similar  way.  We  only  point  out  some  differences*  First,  the  definition  of 
sub-supersolution  is  the  same  except  that  the  boundary  conditions 

BuQ  <  ♦1  <  Bu°  ,  CvQ  <  <  Cv° 

should  be  satisfied  Instead  (1.10),  (1.12). 
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The  main  existence  result  is 

Theorea^l .2.  Suppose  that  Uuq,Vq),  (u°,v°)l  is  a  sub-supersolution 
satisfying  (H^,  (H2).  Then  there  exists  at  least  one  solution  (u,v)  of_ 

-Au  +  o(u)f(v)  a  0  in  Q 
-Av  -  @(u)g(v)  a  0  in  Q 
Bu  -  Cv  -  ♦  on  3Q  , 

such  that  uQ  <  u  <  u°,  vQ  <  v  <  v°.  Moreover,  u,  v  6  W2'p(0)  for  any  p, 

1  <  p  <  4». 

Proof  (sketch).  We  just  give  the  definition  of  the  nonlinear  operator  Tj 
the  other  details  are  very  similar  to  those  for  the  Dirichlet  problem.  For 
(u, v )  6  K,  define  T(u,v)  -  (w,z)  to  be  the  unique  solution  of  the  system 


-Aw 

+  a(w)f (v)  +  w  a  u 

in 

n 

Bw 

=  +1 

on 

an 

-Az 

+  z  =  6°(u)g(v)  +  v 

in 

n 

Cz 

m  lb 
y2 

on 

an  . 

The  existence  and  uniqueness  of  w  and  z  follows  from  Iheorem  II. 1  of  [6] 
(for  z  we  can  also  use  the  results  of  [10]). 

A  result  very  close  to  Lemma  1.3  can  also  he  proved  for  the  boundary 
conditions  (NBC). 

Remark  1.6.  The  operator  -A  in  (NLS)  can  be  replaced  by  two  (possible 
different)  elliptic  second  order  differential  operators  or  even  by  nonlinear 
operators  of  the  form 


■A  u  = 
q 


3- 


with  1  <  q  <  •.  Indeed,  in  this  case  one  can  define  a  nonlinear  operator 
T  by  using  again  (6]  (cf.  also  {15}).  The  more  involved  situation  of  b 
and  c  maximal  monotone  graphs  can  also  be  studied  by  similar  methods. 
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CT1 


§2.  i2iS^2S£2=.2£»2»I*22»»£2£2l* 

In  this  section  we  shall  consider  the  existence  of  a  "dead  core"  for 

solutions  u  of  (NLS),  i.e.,  we  shall  prove  that  the  set 

flQ  *•  {x  e  0  :  u(x)  *  0}  has  a  strictly  positive  Measure  under  adequate 

hypotheses  on  a  and  eventually  on  1*1  _  or  |Q|.  In  fact  much  store 

'  L  OQ) 

precise  information  is  obtained  about  QQ. 

Our  study  will  be  carried  out  by  using  results  concerning  a  single 
nonlinear  equation  but  arguing  in  a  different  way  than  usual  for  the 
combustion  example.  Indeed,  if  (u,v)  is  any  solution  of  (NLS)  (DBC)  [reap. 
(NLS)  (NBC)]  then  u  satisfies 


(2.1) 

-Au  +  f(x)*a(u) 

a  F(x)  in 

n(1> 

(2.2) 

u  *  <P1 

on 

an 

[resp. 

(2.3) 

.  b(u)  -  ♦, 

on 

30]  , 

where  F  S  0 

and  f(x)  «  f ( v( x ) ) 

a.e.  on  SI. 

The  study  of  the  subset  0Q  corresponding  to  solutions  of  (2.1),  (2.2) 
(or  (2.1),  (2.3))  has  occupied  the  attention  of  many  authors,  but,  as  far  as 
we  know,  all  these  results  are  given  for  the  simplest  case  f(x)  =  constant. 
We  recall  the  two  different  approaches  in  the  literature: 

a)  £)  »  RN  ([7])  or  (l  being  an  unbounded  set  ([11]) 

b)  a  being  multivalued  at  the  origin  ([9],  [5],  [20],  [24]). 

More  recently,  a  systematic  study  has  been  made  in  [12]  giving  a  unified  view 
of  both  situations,  but  always  for  f(x)  constant.  Our  results,  in  this 
section,  follow  the  ideas  of  [12]. 


(1) 

Equation  (2.1)  also  appears  in  the  study  of  a  stationary  isothermical 
single  reaction.  (See  [3]  Chap.  3.) 
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2.1  Dirichlet  Problem. 

We  now  prove  parts  i),  ii)  and  iii)  of  Theorem  B.  For  this  we  begin  with 
some  useful  Lemmas. 

a*  2  2  « 

Lemma  2.1.  Let  F  6  L  (ft),  ^  e  C  (30)  and  suppose  that  u  6  H  (0) 

atsxrxrasassiasssx  ■■■■—■  — . m  11 . .  —  111  1 

satisfies 

(2.4)  -Au(x)  +  M2? (x) |u(x) |Psign  u(x)  S  F(x)  in  0 

(2.5)  u  =>*  on  30 

where  f  9  L  (0 ) ,  f  >  0  on  0  and  p  >  0  ^ 1 ^ 0  <  p  <  1  and  0^ 
denotes  the  set 

nx  «  {x  e  n  :  f(x)  >  x)  ,  x  >  o  , 

we  have  the  estimate 


0Q  =  {x  6  0  :  u(x)  “0}  o  {x  6  0^  -  supp  F  s  d(x,3(0x  -  supp  F)  -  (3Q  -  supp 


(2.6) 


Here 
Az  = 


M  m  max{( 
y2X  in  0 


in  -  i 

L  (0)ip 

5  )  , 

Xu2 

,  z  ■  <P  on 


M  „  } 

L  (30) 
30)  for 


for  p  >  0  and  M  =  Izl 

L  (0) 

p  ■  0.  K.  ..  is  given  by 

*  r  H 


(with 


(2.7) 


1 

^  2 1 1 +p+(H-1 ) (p-1 ) ] j p- 1 

Xw2(1-p)2 


Proof.  If  we  denote  by  u+  (resp.  u_)  the  solutions  of  (2.4),  (2.5) 
corresponding  to  the  data  F+,  p*  (resp.  F",  <p~)  then  by  well-known 


(1) 

If  p  ■  0,  equation  (2.4)  should  be  interpreted  in  the  sense  that  there ^ 
exists  w  6  L2(fi)  such  that  w(x)  6  sign(u(x))  a.e.  xeO  and  -Au  ♦  y2fw 
»  P  in  a. 
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comparison  theorems,  we  have  u+  >  0  (reap,  u  <  0)  and  also 
u  (x)  <  u(x)  <  u+(x)  a.e.  x  G  8.  Hence  it  is  clear  that  8Q  3  {x  8  Q  : 

u_(x)  *  0  and  u+(x)  -  0} .  For  the  sake  of  simplicity  we  shall  only  consider 

the  case  F  -  F+,  <p  -  <p*,  the  other  case  being  analogous.  Let  6  h2(Q) 

such  that 

(■  2 

-Au^  +  111  |u^|p  >  F  in  8^ 

(2.8)  \  UX  >  *  °n  30 

u.  >  lul  on  38.  -38  . 

V  A  L  (8) 

We  claim  that  0  <  u(x)  <  u^(x)  a.e.  on  8^.  Indeed1  ,  talcing  F(x)  - 
-Au  +  Xii2up,  it  is  clear  that  F(x)  ■  F(x)  +  Xp2up  -  f{x)|i2u^  and  hence 

<v  2  D  ~ 

F  <  F  on  8^.  Moreover,  -Au  +  Xu  r  ®  F  on  8^  and  thus  by  the  comparison 
results  (cf.  e.g.  [12])  one  has  0  <  u  <  u^.  Therefore  the  conclusion  of  the 
Lemma  will  follow  by  constructing  one  of  such  functions  u^  and  the  set 
{x  6  8^  ;  u^(x)  »  0}  will  give  the  estimate  (2.6)  for  8Q,  we  will  choose 
u^(x)  «  h(|x-xQ|)  for  some  xQ  6  8^.  First,  note  that  for  h  G  C2(E)  and 
any  n  e  (0,1)  we  have 

-Ah  ( |  x-xQ  | )  +  Xu^tlx-XgljP  -  -h"(|x-x0|)  -  (  |  x-xQ  |  )h ' ( I  x~x0 1 )  + 

+  Xu2h(  |x-xQ|  )P  -  -h* ( | x-xQ | )  +  nXii2h(|x-x0|  )P  ♦  (1-n)Xu2h(  |x-xQ|  )p 


(1) 

We  shall  only  prove  this  inequality  for  p  >  0.  If  p  ■  0,  a  natural 
adaptation  of  the  argument  leads  to  the  claim. 


1 


Now,  for  «  fixed  n,  let  be  a  solution  of  the  Cauchy  problem 


(2.9) 


hj(r)  9  hlw2|hn(r)|p»ign  0»n(r>) 

V0)  “  hn<0)  *  0  * 


It  is  easy  to  check  (recall  that  0  <  p  <  1)  that 

2 

(2.10)  h  (r)  -  L  r1-p 

n  n 

where 

(2.11)  Ln  - 


1 

(-.2<14E> 

nin2(i-p)2 


is  a  solution  of  (2.9).  Ne  have 


i£_ 

(1-n)X|i2h  (r)p  -“^h'(r)  -  L  r^K I-iDX^lJ*1 


2(H-1 ) 
1-P 


1  I 


choosing  n  such  that 
(2.12) 
leads  to 


n  < 


-Eli. 


1+P+(N-1 ) ( 1-p) 


-Ahn(|x-x0|)  ♦  Xm\(|x-x0|)P>  0 

for  any  x  9  0^. 

m 

Finally,  consider  the  set  0  ■  -  supp  F.  The  considerations  made 

above  show  that  the  function 

2 

ux(x)  "  Kx,u|x”*o|1’P 

with  K.  given  by  (2.7)  satisfies 
a  fy 

-Au^  +  Xp2uP  >  0  ■  F(x)  in  0 

m 

u^  >  0  -  <p  on  90  n  (IS  ■  supp  <p )  . 


Hence  it  is  sufficient  to  have 

m  m 

(2.13)  u.  >  eax{  <fi,  lul  m  )  on  W  -  (30  n  (30  -  supp  ) ) 

L  (O' 
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to  obtain 


0  <  u(x)  <  u^(x)  on  0  . 

But,  by  the  maximum  principle  we  know  that  u(x)  <  N  on  0  and  this  implies 
that  (2.13)  is  satisfied  if  we  choose  xQ  such  that 

-  llE 

(2.14)  |x-x  |  >  (*“ — )  2 

for  every  x  6  38  -  (3$)  n  (30  -  supp  i)).  The  conclusion  now  follows 
trivially  from  (2.13)  and  (2.14)  (we  recall  that  u^(xQ)  *  0). 

I 

Statements  ii)  and  iii)  of  Theorem  B  follow  immediately  from  the  above 
Lemma.  We  remark  that  the  constant  K.  given  in  (2.7)  is  such  that 

A  f  y 

K.  0  when  X  ^  0  or  0  and  that  K.  if  X  7*  or 

A  f p  *»P 

U  ?  +*.  Then,  if  0^  is  bounded  and  not  empty,  estimate  (2.6)  shows  that  the 
measure  of  8Q  is  positive  at  least  if 

•  !iE 

diameter  of  (0^  -  supp  F)  >  (“ — )  2 

X,u 

(assumed  0  4  p  <  1).  Therefore,  if  8^  is  given,  0Q  "exists"  if  M  is 
large  enough  or  M  is  sufficiently  small.  In  the  simple  case  of  problem 
(0.1),  (0.3)  with  f(s)  -  sp,  0  <  p  <  1,  it  is  easy  to  find  a  critical  value 

of  m  (now  depending  on  p,  Y  and  diameter  of  8)  such  that  8Q  is  not 

empty  if  u  >  Uc.  When  H  -  1  direct  computations  show  that,  for  p,  Y  and 
8  fixed,  the  function  u  is  strictly  positive  if  w  <  Uc  (see  e.g.  [3]  and 
(221). 

we  shall  prove  part  i)  of  Theorem  B.  Indeed,  we  shall  prove  that  if 
p  >  1  then  for  any  value  of  X,  w  and  diameter  of  8  thar*  exist  functions 

(u,v)  satisfying  (DP)  (with  a(s)  ■  |m|^'m)  such  that  u(x)  >  0  on  8.  To 

do  this  we  shall  consider  the  worst  of  the  cases  i.e.  when  diameter  (8)  »  +* 
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(for  instance  N  *  1  and  (2  =  (0,*))  and  even  for  a  larger  class  of 


I 

I 


! 

I 

! 

i 

i 


nonli near! ties  a. 

2 

Lemma  2.2.  Let  u  S  H  (0,“)  satisfying 


(2.15) 


-u"(x)  +  f(x)a(u(x))  a  0  ,  x  6  (0,«) 

u(0)  »  1  , 


where  a  is  a  maximal  monotone  graph  snch  that  0  6  a(0)  and  the  function 
j(s)  *  /g  a°(r)dr  satisfies 


(2.16) 


'  j  (a ) 


(These  hypotheses  are  satisfied  when  o(s)  *  | s | p  >  i.)  Assume  that 
"  •  " 

f  6  L  (0,®)  and  0  <  f(x)  <  m 2  a.e.  x  6  (0,*),  for  some  Then 

u(x)  >  0  for  any  x  6  [0,®). 

Proof.  We  shall  use  some  ideas  of  (7)  and  [15].  By  reasoning  as  in  the  proof 


of  Lemma  2.1  we  can  always  suppose  without  loss  of  generality  that 


-1 


a  (0)  »  0  and  that  a  is  singlevalued.  By  a  comparison  argument  completely 

.2 , 


analogous  to  the  ones  in  the  proof  of  Lemma  2.1  we  show  that  if  u  8  H  (0,®) 
satisfies 


-u"(x)  +  m2a(u(x))  *  0  on  (0,®) 
u(0)  -  1 


then  u(x)  <  u(x)  for  any  x  8  (0,+“).  Thus  it  suffices  to  prove  that 
u(x)  >  0  for  any  x  8  (0,®).  Suppose  that  u  has  compact  support  and  we 
shall  obtain  a  contradiction.  The  maximum  principle  implies  0  <  u(x)  <  1 
and  hence  u"  8  L  (0,®).  Thus  u  8  cVtO,®))  with  u*  >  0.  Let 
R  ■  sup{x  :  u(x)  0}  (R  >  0  and  finite  by  assumption).  As  u'(R)  *  0  it  is 
not  difficult  to  see  that  u'(x)  <  0  and  u(x)  >0  on  (0,R)  (it  is  a 


consequence  of  u"  >  0).  But  (2.16)  yields 


u(R) 


ds 


'T7F) 


rR  H'<*> 

*/  o - dr  -  +- 


/j(u(r>) 


and  we  will  get  a  contradiction  by  estimating 


w(x)  m  (u*(x))  we  have 


u'(r) 

/j(u(r)) 


on  (0,R).  Defining 


(j(u))*  -  o(u)u‘  =  1-  u*  u"  -  l-~-  (u')2r  . 

-  -  -  “2  ”  -  “*2  ~ 

But  w(R)  -  0  and  j(u(R))  ■  0.  By  integrating  we  get  j(u) 


w,  and, 


finally. 


a  contradiction. 


u'(r)  _ 

/*  - - dr  -  /ST  /*  d.  <  «.  , 

°/jM77>  20 


Remark  2.1.  By  arguing  in  a  similar  way  as  in  [11]  we  can  prove  ([14])  that 
if  U  is  an  unbounded  subset  of  R1*,  the  maximal  monotone  graph  a 
satisfies 


(2.17) 


and  u  satisfies 


J o  ;-d3  <  (j(»)  -  /«  <*°(r)dr)  , 

✓j(s)  0 

-Au  +  f(x)o(u)  8  F  in  8  (f  »  A) 
u  »  i p  on  30 


where  F  and  v  are  assumed  with  compact  support,  then  u  has  compact 

support.  We  point  out  that  the  improper  integral  (2.17)  converges  when 

o(s)  -  |s|p  sign  s  if  and  only  if  0  <  p  <  1  and  hence  the  compactness  of 

the  support  of  u  is  an  obvious  consequence  of  Lemma  2.1.  On  the  other  hand, 

we  do  not  know  (for  N  >  1 )  if  under  the  general  hypothesis  (2.17)  the  set 

0  has  positive  measure  for  0  bounded  and  perhaps  I  FI  m  and 

L  (0) 
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sufficiently  sea 11  depending  on  the  size  of  0 


M  . 

L  OQ) 

Remark  2.2.  Lemma  2.1  (and  then  Theoree  B)  can  also  be  obtained  when  the 

operator  -A  in  (NLS)  is  replaced  by  other  elliptic  second  order  differential 

operators  as  in  Remark  1.6.  The  new  definition  of  the  functions 

u^(x)  «*  h(|x-xQ|)  in  the  proof  of  Lemma  2.1  can  be  found  by  the  methods  of 

[12J. 


2.2.  Nonlinear  boundary  conditions. 

Statement  iv)  of  Theorem  B  will  follow  as  in  the  preceding  section  by 
considering  the  nonlinear  equation 


(2.18) 


(  -Au  +  |i2f(x)|u|^  sign  u  9  P  in  0 
Bu  ■  +  b(u)  ■  ♦  on  30 


where  f  €  L  (0),  f  >  0,  0  <  p  <  1(  b  is  C1  non-decreasing  with  b(0)  -  0, 
FeL*(0)  and  ♦  8  C2(3fl). 

First,  we  remark  that  "interior  estimates"  for  0Q  can  be  obtained  as  in 
Lemma  2.1.  More  precisely,  we  have 

*  JbE 

(2.19)  Qq  a  {x  6  -  supp  F  s  d(x,3(Qx  -  supp  F) )  >  — )  2  } 

*  (1 ) 

where  now  M  *  lul  .To  show  this  it  is  sufficient  to  choose  x0  in 

L  (Q) 

•  “ 

such  a  way  that  u^  >  M  on  30,  being  0  ■  0^  -  supp  F  in  the  proof  of 
Lemma  2.1.  It  is  clear  that  (2.19)  does  not  give  any  information  about  the 
behavior  of  0Q  near  the  boundary  of  0^  -  supp  F. 


One  obtains  estimates  for  M* 
Lemma  3  in  (12)). 


by  means  of  comparison  theorems  (see  e.g. 


To  improve  the  estimate  (2.19),  we  introduce  the  following  notation: 
given  a  smooth  curve  r  in  and  xQ  6  we  define 

(2.20)  0(xQ,r)  *  inf{cos(n(x),  x-xQ)  :  x  6  T} 

where  n(x)  =  (n^ (x),. . . ,n  (x) )  is  the  unitary  outward  normal  vector  to  T 

at  x  and  (n(x),  x-xQ)  denotes  the  angle  between  the  vectors  n(x)  and 

x-xQ.  It  is  clear  that  the  value  of  0(xQ,r)  depends  essentially  on  the 

"geometry"  of  T.  If  for  instance  r  =  30  and  (1  is  a  convex  bounded  set 

of  RN  it  is  easy  to  see  that  0(x^,r)  >  0  when  x^  6  0. 

2  00 

Lemma^^S.  Assume  that  u  6  H  (0)  n  L  (0)  satisfies  (2.18).  For  X  >  0, 

let  0^  =  {x  6  $2  :  f(x)  >  X}.  Moreover,  suppose  0  <  p  <  1  and 

(2.21)  0(xQ,  ®^x  ~  SUPP  n  3(2)  >0  VxQ  6  0^  -  supp  F  . 

Define 

r  m  3  (0^  -  supp  F)  n  3(2  n  supp  ^ 

Then 

0Q  a  {x  6  0^  -  supp  F  :  d(x,F)  > 

d(x,  3 (£2^  -  supp  F)  -  3 £2)  >  — )  2  } 

X,u 

* 

where  M  =  I ul  _  • 

L  (£2) 

P-oof .  Arguing  as  in  Lemma  2.1  we  only  consider  the  case  F  >  0  and  ♦  >  0. 
Let  0=0^-  supp  F.  By  using  again  comparison  results  (cf.  e.g.  [12))  it  is 
not  difficult  to  see  that  if  u^  satisfies 


(2.23) 

-Au^  +  Xy2uj^  *  0 

in 

£2 

(2.24) 

UX  >M* 

on 

3£2 

-  30 

3u, 

(2.25) 

3n  L  (3£2) 

on 

r  - 

30  n  30  n  supp 

3uX 

(2.26) 

x— —  >  0 

3n 

on 

3  £2 

n  (30  -  supp  ♦) 

(1-p)»*l 


L  ( 3£2) 


2Kx,u0(Vr) 
*  JL-.E 


izE 

i+p 


and 
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then  0  <  u(x)  <  u^(x)  for  x  9  ».  From  the  proof  of  Lemma  2.1  we  know  that 


the  function 


satisfies 


Vx)  -  kaJ*-xoI 


-Au.  +  Xy  u?  >  on  Q 


for  any  xQ  9  Q.  Condition  (2.24)  is  satisfied  if 


(2.27) 


|x-x  I  >  — )  2  Vx  9  30  -  30  . 

x,u 


On  the  other  hand, 
3u,  N  3u 


a^-(x)  “  a7^(x)*ni(x)  "  ka<w(T^)  Ix-x0I  "P  *»(*(*>'*-*0>  > 


i±E 

>  |x‘xol1_P  0(xo'3“  n  3Q)  • 


Thus  (2.26)  is  a  consequence  of  (2.21),  and  (2.25)  holds  if  we  choose  xQ  6 


satisfying 


(1-p)l*l  1-p 

|x.x  |  „  f _ OQ)'1+P  r 

lx  V  >  l  2K  0{x  ,T)  >  Vx  9  r  . 


This  completes  the  proof. 


Remark  2.3.  Part  iv)  of  Theorem  B  follows  from  Lemma  2.3  if  we  set  F  =  0» 
(2.21)  holds  easily  if,  for  instance,  f(s)  >  >  0  Vs  6  R  (as  in  the 

combustion  system)  and  (1  is  a  convex  set. 

Added  in  proof:  After  the  completion  of  this  work,  the  authors  have  learned 
that  C.  Bandle  and  I.  Stakgold  have  obtained  recently,  results  similar  to  our 
Theorem  B  for  the  particular  system  (0.1),  (0.3),  by  using  different 
methods.  On  the  other  hand,  some  results  related  with  Remark  2.1  can  be  found 
in  a  forthcoming  work  of  M.  Schatzman  ("Stationary  solutions  and  asymptotic 


behaviour  of  a  quasi linear  degenerate  parabolic  equation"). 


-2 
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